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Abstract 

In this paper we describe the theory of a device that is able to extract 
energy from an external magnetic field. The device is a cylindrical mag- 
netic insulator that once put in rotation makes electromagnetic angular 
momentum to be stored in the electromagnetic field in contrary direction 
to the mechanical angular momentum of the device. As a consequence due 
to total angular momentum conservation the device increases its angular 
velocity (when e/i > 1 j3 and all conservation laws are rigorously satisfied. 
The voltage generated by the device is found solving explicitly Maxwell 
equations for rotating magnetic insulators in external fields a subject that 
have provoked lots of polemics in the literature and which we hope to be 
here clarified due to our pedagogical presentation. 



1 Introduction 

In this paper we introduce the theory of a device that can in principle ex- 
tract energy from an external magnetic field. The idea for the device came to 
us from the theoretical analysis of the 1913 Wilson & Wilson experiment 
(WWEjl described in Section 5. Since the correct theoretical explanation of 
the experiment gave rise to controversies in the literature due to difficulties 
in properly applying Maxwell theory in rotating framed we decided to give a 



•email: capelas@ime.unicamp.br 
^Natural units are used in the paper. 

■^The experiment has been repeated by Hertzberg and collaborators [5] in 2001. 

sample may be found in 1111 1131 1171 1181 120| and the inumerous references to related 
issues in those papers. 

■'See, e.g., the most quoted SchifT's paper 1191 where it is (wrongly) claimed that one need 
use General Relativity to explain electrodynamics phenomena in rotating frames. Besides that 
let us recall the almost continuous difficulty that some have in each generation to understand 
simple electromagnetic phenomena when moving boundaries are present, as e.g., in the appli- 
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very pedagogical presentation for the problem of expressing Maxwell equations 
in a coordinate invariant way using the theory of differential forms. Then, using 
that formalism and the correct jumping conditions for the electromagnetic field 
variables between moving boundaries separating two different media we get in 
a very clear way the electromagnetic field as measured in the laboratory in the 
wwe| We think that our approach clears up immediately what is wrong with 
some other attempts to find a correct description of the WWE in the sample 
of papers quoted. The description of the machine to extract energy from an 
external magnetic field is given in section 6 where it is observed that the ma- 
chine once puts in motion makes electromagnetic angular to be stocked in the 
electric plus magnetic fields in the opposite direction of the mechanical angular 
mechanical of the device (when > f) and thus due to angular momentum 
conservation the machine has its mechanical angular momentum increased. Of 
course, all the conservation laws are in operation in our device and the energy 
being generated comes from the external electromagnetic field. 

2 Maxwell Equations 

The formulation of Maxwell equations in intrinsic form requires as preliminary 
mathematical structures: 

(i) a 4-dimensional orientablcd manifold M, the bundle of non homogeneous 
differential forms f\T*M = ®lZtKT*M, where /\''T*M is the bundle of r- 
forms and 

(ii) the differential operator d : f\''T*M ^ /\''^^T*M, (P = 0. 

Indeed, Maxwell equations deals with a field F E sec /\^T* M which is exact, 
i.e., F = dA, where A e sec /\^T*M and a current J e sec /\^T*M that is also 
exact, i.e., J = dH, where H e sec /\^T* M . The set of equations 

dF = 0, dH = -J, (1) 

is known as Maxwell equations and of course, we have 

dJ = 0, (2) 

i.e., the current is conserved. 

Maxwell equations are invariant under diffeomorphisms as it happens with 
all theories formulated with differential forms. This means that if ft, : A/ — > M 
is a diffeomorphism, then denoting as usual the pullback mapping by h*, since 

cation of Faraday's law of induction in the homopolar generator. A thoughtful discussion of 
the application of Faraday's law in its many equivalent formulations may be found in |16) . 

^Application of differential forms formalism to the WWE already appeared in [l][6]. Our 
approach details each step in the derivation and hopefuly is a little bit more pedagogical. 

°By orientable we mean that there exists in M a global 4-form Tg G sec /\'^T*M. If the 
manifold is not orientable then it is necessary to use besides the concepts of pair form fields 
also the concepts of impair form fields. Details may be found in [6]. See also a thoughtful 
discussion in |14) . 
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the differential d commutes with the pullback, i.e., dh* = h*d we have that the 
fields F' = h*F,H' = h*H,3' = h*3 satisfy 

dF' = 0, dH' = - J'. (3) 

Note that we did not make until now any requirements concerning the topol- 
ogy of manifold M, so to proceed we take arbitrary coordinates {x^) covering 
U C M and denoting (e^ ~ djdx^ = dp) the corresponding coordinate basis 
for TU and by (i?'^ — dx^) the basis for T*l] dual to the basis (e^) we write 

F =\^Fp,'d^ N-d\ 

H^^Hp^^^Ar, (4) 

Maxwell equations are supposed to describe the behavior of electromagnetic 
fields in vacuum or in material media. But whereas, according to Feynman, 
the field F is fundamental (it define the Lorentz force acting on probe charged 
particles moving in the field) the field H is phenomenologicajl. In fact, if we 
write Fourier representation for F and iJ, 

Fp,{x) = /d4A;F^,(fc)e-^('=o-"+'=i-^+'=^-'+fc3-^) Jd^k^^,{k,x), 
Hf^^ix) = Jd^k ^^^(fc)e-»('="-''+'^i-'+'=^-'+'=^-') Jd^kSjp^ik, x), (5) 

in general we have, defining the 2-form valued function ><r : Mx sec 
sec/\ T* M (the frequency constituent equations of the medium) that 

^{kf),x) = >c{^{ko,x)). (6) 

The function x may even be in some media a nonlinear function of ^ (see, 
e.g., [2]) but in what follows we will consider only non dispersive media {x is 
independent of fcg) in which case we can define a linear constituent function (an 
extensor field) 

X : sec /\^T*M sec /\^T*M, 

H = x{F). (7) 

In this case we have for the components H^^ of H , 

H,. = y^f.Fap. (8) 

^Take notice that coordinate vector fields are denoted using a bold symbol d, e.g., d/dxf' = 
dfj_, whereas we use the symbol d/dxi^ := to denote the usual partial derivatives. This 
means the following: Let / : M — >■ R a differentiable function and let {U, be a chart of the 
atlas of M such that for e G (7, ij{e) = (a;°, a;^, x^, i^). Then, qItt/ = gf^r/ o ip. 

*At least at the classical level. This is so because the calculation of H needs in the general 
case sophisticated use of quantum theory. 
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Of course, we have the obvious symmetries for the components of the con- 
stituent extensoi@ >c, 



2.1 Enter (M,^,V,Tg,t) 

We can proceed with the formulation of Maxweh theory without the introduc- 
tion of additional mathematical objects in the manifold M . If you are interested 
in knowing the details, please consult [5]. From now on we will suppose that 
electromagnetic fields are to be described in a Lorentzian spacetime structure 
(M, g, VjXgjt), where {M,g) is a Lorentzian manifold, V is the Levi-Civita 
connection of g, Tg = ^J\Aetg\if' A A i?^ A is the metrical volume element 
and t denotes that the structure (Af, g) is time orientable. 

Remark 1 As it is well known a spacetime structure (M, g^V ,Tg,'\') when the 
Riemann curvature tensor of V is nan null represents a gravitational field gen- 
erated by an energy-momentum tensor T in Einstein's General Relativity ( GR) 
where Einstein equation is satisfied. R is a simple exercise to find an effec- 
tive spacetime structure, say (Af, g, V, Tg, 'f), to describe some material media. 
However, such .structure has nothing to do with GR. 



2.1.1 Reference Prames 

We shall need in what follows the concepts of reference frames, observers and 
naturally adapted coordinate systems to a reference frame in a general Lorentzian 
spacetime structure (Af, g, V, Xg, f)- 

We define a reference frame in U C Af as a unit timelike vector field Z 
pointing to the future. We have 

g(Z,Z) = l. (10) 

An observer is defined as a timelike curve a in M (parametrized with proper 
time) and pointing to the future. We denote by a^s the tangent vector at cr(s) 
and 

g(crH.s,cr*s) = 1- (11) 

We immediately realize that each one of the integral lines of Z, is an observer. 

Finally we say that coordinates (cc^) covering U G M is a naturally adapted 
coordinate system to a reference frame Z (nacs|.Z) if 

A detailed classification of reference frames in a Lorentzian spacetime may be 
found, e.g., in |15j . For a classification of reference frames in a Riemann-Cartan 
spacetime, see, e.g., [3]. 



^We recall that these are the same symmetries of the components of the Riemann tensor 
of a metric compatible connection. 
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2.1.2 Formulation of Vacuum Maxwell Equations in (Af, g, V, Tg, t) 
Given the structure {M,g,V ,Tg,^) it is convenient to define 

G = ir-^H (13) 
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where ★ : sec /\^T*M sec/\ ^T*M is the Hodge star operator and -k ^-k = 

9 9 9 

irk^^ = Id. It is also a standard practice to introduce J £ sec /\^T* M such that 

99 

J = * J. (14) 
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X ■■ sec /\^T*M sec /\^T*M, (15) 



We define also the (ex) tensor field 

I, 

Writin£^ 

G'=^G^.i?''Ar = iG'"^i?^A^, (16) 

we have 

G,. = lx;fFo.p, G^'' = \x^'''^^F^P (17) 
and of course, the components of the constituent extensor x satisfy 

In vacuum the relation between the fields G and F is 

G^F (19) 

and we write Maxwell equations as 

dF = Q, d*F=-*J. (20) 
9 g 

or better yet, applying the inverse of the Hodge star operator to both members 
of the non homogenous equation we can write 



dF = 0, 5F ^ -J 

9 



We recall thalR 



= ^ -v/kietg|g^i''i...g'^'''^''e^, ^ r^-^^-"- . (21) 

9 (n - p) \ 



■"-"(j?^) is the reciprocal basis of {"Q^), i.e., = S'i . 
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Then, write 



H = ^F= ^Hp^r ^ ^HP'^^p A (22) 



and let us calculate -kF. We have 

9 



1 



9 2 g 
1 1 



= --F^.V^drt^g^^g'^'^e^^p,,?'' A i?^ (23) 



From Eq.(122]) and Eq.(|23]) we have 

Hp„ = iv/|detg|5''"ff"^ea^p.F^. (24) 



2 

or taking into account that 



= -V^detg|e;r^^M^, 



i^p. - *Gp, - -^\detg\ep^!'.Gp, (25) 

we get 

= gP^g'^'-Fp, = ^{g^^g'''' - g'^'g^nF^.u- (26) 

Comparing this equation with Eq. p7|) gives for the components j^p"'^'^ of 
the constitutive (ex)tensor of the vacuum 

;tP-M- ^ (gPMg- _ gP'^g^f^y (27) 

Remark 2 Before proceeding let us recall that the similarity between the com- 
ponents of ^P'^t^'^ and the components Pf^P-^ of a Riemann curvature tensor R 
with a constant Riemann curvature scalar K . Indeed, we have 

R""'''' = ^if^f'' -f'g^n (28) 

and thus K = R = SpR".a^. Take notice that R is not the Riemann curvature 
tensor of the Levi-Civita connection Dofg. 

Exercise 3 Show that in components Maxwell equations dF — 0. dH = — -k J 

9 

when K has its vacuum value read 

dpFap + dpFp^ + da^Fpp = 0, (29) 
^ =ap(v^d^F^^) = r. (30) 
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Solution 4 We have 

dF = ^d{Fa0 A^" A^^) 

= ^{df^F^p + dfiF^,^ + do.Ffi^)^^' Ad'^Ad^ 

and Eg. ([29)) follows. 

Next we define J — J^d'^ := pi?" — jid"^ and get 

★ J = ^ V^d^JA5^"eKpa/3i9^ Ad^Ad^. (31) 
Then the equation dH = — -k J reads 

9 

^{d^H,,p + dpH^^ + d,,Hpf,)d^ A I?" A i?/^ = -ividetg|J'^e«^„^i?^ A i^" A 

(32) 

and 

{df,Hap + dpHfj,a + daHfSf,) = - V|dctg| J^eKAiQ/j. (33) 
Multiplying both members of the last equation by e^'^"^ and recallinS^ that 

e^"°''eKpa/3 = 3!5;J (34) 



we 



have 



e^--^{d^H^p + + 5„F^^) = -3!V|dctg|j\ (35) 

= 5,,(/^J^r^'^^) 
= 29^ (/^^F^^) 

= -29^(/^^F^^). (36) 

Analogous calculations give that the other two terms in the first member of 
Eq.dSil) are also equal to -2^^{^J\Aetg\F^'^). Then Eq.dSS]) becomes 

-y^d,{./W9\F^^^) = J' 

V|detg| 



and Eq.((30|) is proved. 

^^See, e.g. page 111 of [8]. 
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Remark 5 We present yet another simple proof of Eg. pop which however pre- 
sumes the knowledge of the Clifford calculus. We introduce the Dirac operator 
d — '^^'Va/^x^^ and recall that dM = d A M for any M e sec /\T*M ^ 
secC£{AI,g) where C£{AI, g) is the Clifford bundle of differential form^^. We 
have 

dH^d*G = dA*G = -*J 

9 9 9 

Now, we calculate d A-kG = d A (GTg) in orthonormal hasisV^Ba^ for TU and 

9 

dual basis {6''} for T*U with U C M. We have 
d*G^dA{GTg) 

= -0'^A(VeJGr<,)) 

= (0'^(Ve,G)xg + rG(Ve,rg) + (Ve.G)rg0'^ + G{V,^Tg)e^) 
= -\ (^"(Ve.G) - (Ve.G)r)rg - i( - e'^GiVe^Tg) - G{Ve^Tg)e-^ 
But 

-- {e'^iv^G) {v^^G)e-)Tg) = -{d.G)Tg 

and 

Ve.Tg = Ve,(0° AO^AO^A e^) = -r°,o^a - r^^iTg - r^^^r^ - T%Tg = o, 



because 
Then, 



or 



Now, 



r°aO — r^al — ri'al — ^^al — 0- 
-{djG)Tg = ~jTg 

BaG = J. (37) 



djG = ii?''j(I?^G"'^i?„ A - D^G^'^'&p 



{d,G^^ + V^.^G'^P)^^ = {d.G^P + ^^^^G^f')d, 



^The ClifFord product of Clifford fields is denoted by justaposition. 



"Recall that e'^Cet) = &f. Also g(ea, et) = »)ab and g(0'*, e'') = ry^^ = ■ 0^ ^ fla^^b 



8 



and we get recalling that in vacuum G = F 



which is Eg. pop again. 



3 Maxwell Equations in Minkowski Spacetime 
3.1 Vacuum Case 

We next suppose that electromagnetic phenomena take place in a non dis- 
persive material medium in Minkowski spacetime, i.e., the structure M = 
(M'*, g, V, Tg, t). Since the Riemann curvature of the Levi-Civita connection V 
is null (i.e., fl(V) — 0) there exists global coordinate^ (x'') such that denoting 
by (e^ = ^/^yi^') a basis for TM and (7^ = dx'^) the basis of T*M dual to (e^) 
we have 

g = ® 1" (38) 

where the matrix with entries 77^1/ being the diagonal matrix diag(l, —1, —1, —1). 

In this case the components of constitutive (extensor) x of vacuum in 
the (nacs|/) {x'^} to the inertial frame L = = d/dx'^ are 



Thus 



(39) 



(40) 



Denoting by (F^,,) and (C" 
QP'^ we have 





El 


E2 


E3 \ 


-El 





-B3 


B2 


-E2 


B3 





-Bi 


V -E3 


-B2 


Bi 


/ 



respectively, the matrices with elements F^i^ and 
/ -Bi -B2 



( ' 


-Bi 


-B2 


-B3 


\ 


Bi 





E3 


-E2 




B2 


-E3 





El 




V B3 


E2 


-El 








Of course, Maxwell equations in coordinates in the Einstein-Lorentz Poincare 
gauge reads 



^ 

dxP 



and 



dxP 

dF''" 
dxP 



5x° 

r. 



-0, 



(41) 
(42) 



^^These coordinates are said to be in Einstein-Lorentz- Poincare gauge. Note that (x**) is a 
(ncsa|eo) where bq = L (the laboratory frame) is an inertial reference system, this adjective 
meaning that Veo = 0. 
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3.2 Non Dispersive Homogeneous and Isotropic Lineeir 
Medium Case 

We suppose in what follows that a linear non dispersive homogeneous and 
isotropic medium (NDHILM) is at rest in a given reference frame V inU C M 
which has an arbitrary motion relative to the laboratory frame that is here mod- 
elled by the inertial frame L = Bq = d/dyp. Let (x''^) be coordinate functions 
covering U that are (nacsjy) and {d/dx'^^) be a basis for TU and (i?''' = dx'!^) 
the corresponding dual basis for T*M. Writing in this case 



9 = 9'^,^'"'^^"' = Vi..rAl^'' 

we have 

V = v'l'd/dx'^ = -^d/dx'° = va/ax" = v° 



^500 



d 



r 



d 



(43) 



(44) 



We now write 



G 



1 



jGpa7^A7'^ = V%A7<,. 



(45) 



We use the following notations for the elements of the matrices {Fpa), {Fpa) 
[Gp^) and [G'p^) as 



(F.a) = 



(G,a) 








El 


E2 


E3 


\ 




-El 





-B3 


B2 




"E2 


B3 





-Bi 




V 


-E3 


-B2 


Bi 





/ 


/ 





Di 


D2 


D3 


\ 




-Di 





-H3 


H2 




-D2 


H3 





-Hi 




\ 


-D3 


-H2 


Hi 








(G' ) 



/ 





E[ 


E'2 


E's \ 




-E[ 







B'2 











-B[ 


V 






B', 


/ 










(46) 


/ 





D[ 


D'2 


D', \ 




-D'l 




















V 






H[ 


/ 



(47) 



Then Maxwell equations read in the coordinates (x''^) (which are a (nacsjV)): 

dFL, dFl^ dF^^ 
ax'" dx'P dx'" 



and 
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3.3 Minkowski Relations 

We want now to define using differential forms the concepts of electric and 
magnetic fields and induction fields in a NDHILM. 

Given an arbitrary reference frame Z in Minkowski spacetime let z — g(Z ^ ) 
be the physically equivalent 1-form fielcQ. 

Definition 6 The electric E and the magnetic B l-form fields and the D and 

z z z 

H 1-form fields according to the observers at rest in Z ara^T 



E := zjF, B := zj-k F, (48) 

Z Z 

D=zjG, H:=zjirG. (49) 

z z 

We immediately have 

F = zAE-*{zAB), (50) 

Z Z 

G ^ zAD ~i.{z AH). (51) 

z z 

Let V be an arbitrary reference frame in Minkowski spacetime {v — giV ,)) 
where the NDHILM is at rest. 

Definition 7 The electric D and the magnetic H 1-form fields according to 

V V 

the observers at rest in Z are related with the E and B by: 

V V 

D eE and H -B. (52) 

V V V ^ V 



We immediately have 



G = - (e/x - 1) [vA{v^F)] + -F. (53) 

which will be called Minkowski constitutive relation. 
Exercise 8 Prove Eg. (15^. 

Solution 9 From Eq. ([STjl and Definition [7] we have 

G = ev A{vjF)- -ir[v A{vjirF)] (54) 

Now using the identities Eq. (|137|) and Eqs. (jl38p in the Appendix we have 
*[v A 7^ F)] i^[v A A F)] = ★ * [v_i{v A F)] = -v_i{v A F) 



-{v_iv)F + v A {v_iF) 



^^We will also call z a reference frame. 

^^The symbol * denotes the Hodge star operator in Minkowski spacetime. 
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and then 

G^evA {vjF) + -[F-vA {vjF)] 
M 

= -{en~ 1) [vA{vjF)] + -F. 
M M 

3.3.1 Coordinate Expression for Minkowski Constitutive Relations 

We want now to express Minkowski constitutive relations in arbitrary {x^} 
coordinates covering U C M and let Fx^j,, G^^ and be the components of F, 
G and the reference frame V in an arbitrary natural bases {c?^ = d / dx^} and 
{•d^ = dx'^}. Let us calculate vaG. We have 

vaG = ie --\ va[vA[vaF)] + -vaF. 
V M/ A* 

From identity Ea. (|137p in Appendix we have 

vj[v/\{vjF)] = {vjv){vjF)] - t> a a {vjF)) 

= VjF ~ {v /\v) A [vjF)) 

= vjF. (55) 

Then 

vjG = evjF (56) 



and in components Eq. (|56p is 

and since (iCa^jz?" A z?'^) = G^^^d", i^^j (i^^^i?" A iS^^) = i^^^i?" we have 

v^G,,, = ev^F^,. (57) 
Let us now calculate v_i ★ G. First we note that 

★G= fe--)* bA(T;jF)l + - ★ 
and using the same steps as the ones leading to Eg. (1551) we first get 
v_i * [vA{v_iF)] = ^vA{v A (vjF))] = A v)A{v_iF) = 

and then 

vj*G=-vj-kF. (58) 
A* 
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Now, 



2 V2 
and 

and finally Eg. dSSp reads in components 

fiv-C^^e^,^^ = v-Ff^^e^,^, (59) 

or in equivalent forms 



3.3.2 Polarization, Magnetization, Bound Current and Bound Charge 
Fields 

We define moreover the polarization 2— form field 11 by 

n-.F-G. (61) 
Given an arbitrary frame Z we decompose 11 as 

Il:^zAP-*{zAM), (62) 

z z 

where P and M are called respectively the polarization 1-form field and mag- 

z z 

netization 1-forni field in Z. Moreover, from Ea. (l52t and Ea. (l53l) we get 
P = E D and M = B H. 

z z z z z z 

From the non homogeneous Maxwell equation d-k G — — J we have 

d*F = -J + d*n. (63) 

The field 

J:^-d*U (64) 

is called the bound current 3-form. Given an arbitrary frame Z we decompose 
J as 

J:^zAj + p (65) 
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where J is called the bound current 2-forni field and p is called the bound charge 

z z 

3-forin field. 

In the coordinates (x'') and (x'^) respectively adapted to the inertial labo- 
ratory frame L = Bq and to an arbitrary frame V we write 



(66) 



where the entries of the matrices (lipa) and (Hpo.) are, respectively, 



/ 
Pi 

P2 
\ P3 



-Pi 


Ms 
-Ma 



-P2 
-Ma 


Ml 



-Ps 

Ma 
-Ml 





(Hp.) := 



/ -Pi -P^ 
P[ -M^ 



-PL \ 

J 
(67) 



3.4 Constitutive Relations in a Uniformly Rotating Frame 
V. 

We introduce besides the coordinates {x''} also cylindrical coordinates (x^-) in 
M.^ — {0} which are also (nacs|L = cq). As usual we write 

xO = x°, 

x^ = cosx^, 

2 1-2 

x = X Sin a; , 



X 



3 3 



a;^ (68) 



We will simplify the notation when convenient by writing {x°, x^, x^, x"^} := 
{t,x,y,z} and {x^ , , x"^ , x^} := {t,r,(l),z}. 

Next we introduce a particular rotating reference frame V € sec TM in 
Minkowski spacetime: 

where 

V := ^v2 + v2 = wr, (70) 

with Lo the (classical) angular velocity of V relative to the inertial laboratory 
frame L = e^. 
Since 



14 



we can also write taking into account that 7 = 1/ Vl — v^, 
^ Id oj d 

= 'y(^dt+v^d^^ (72) 

and 

V := ) = 7 (rfi + vr#) . (73) 

Remark 10 It is important to recall that although the coordinates {t,r,(l),z} 
cover R** — {0} the reference frame V if realized by a physical system can only 
have material support for r < 1/w. 

We introduce next the vector field 

V = — wy^x + wx9y 

= v^a^ = ve^ (74) 

It is quite obvious that the vector field v represents the classical 3- velocity 
of a (material) point whose 3-dimensional trajectory in the spacelike section 
determined hj L — Bq has parametric equations {r cos ut,r sin cot, zq) where zq 
is an arbitrary real constant. 

Take notice that in engineering notation the vector field v is usually denoted 
by w = ve^. Wc will use engineering notation when convenient in some of the 
formulas that follows for pedagogical reasons. 



A (nacsjV) Before we proceed wc introduce explicitly the transformation law 
between the coordinates (x^) that are (nacs|eo) and {x''^) that are (nacs|y). 
These are 



xO = x'O, 

x^ = x'^ coswx^ — x'^ sinwx^, 

x^ = x'^ sinwx^ + x'^ cosojx^. 
,3 /3 



x^ = x'^. (75) 

To simplify the notation we will write {x'^, x'^,x'^, x'^) := {t' , x' , y' , z'). 
We define also cylindrical coordinates (f*) naturally adapted to V by 



x'° = (76) 

x'^ = x'^ cosf'^, 
x'^ = x' ^ sinx'^, 
x'^ = x'\ 
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We moreover simply the notation writing a;'^, x'^, x'^) := (i', r' , (j>', z'). 

Now, it is trivial to see that 

t' ^t, r' ^r, (j)' ^ (jj- ujt, z' = z. (77) 

This has as consequence the obvious relations 

d£ = dt, dr' = dr, dtf)' = d(t> - udt, dz' = dz (78) 

and the not so obvious ones 

d d d d d d d d d 

dF^di^'^d^' d?^d^' d^^d^' 'd?^d^' ^ ' 

Then we see that 

VI - v2 at' ^goo dt' 

which shows that indeed (x'^) is a (nacsjZ'). 

We also notice that since r' — r and ~ '§4, '^^'^ write Eg. ([7^ aj^ 

Id Id , d , d 

r dcp r d<p dx dy 

9 ,9 



~ ^dx'^" ^dx'^- 

We write the metric g using the coordinates (a;'') and {x^) as 



g = g'f.j.dx''^ «) dx"" 

= (1 - uj'^r''^)dt' (g) dt' + 2uy'dx' (g) dt' - 2ux'dy (g) dt' - dx' (g) dx' - dy' (g) dy' 
-dz'(g)dz', (82) 

or 

g = gf^^dx'^" (g> dx"" 
= (1 - v^)dt g)dt- 2ujd(j)' g)dt- dr' <g) dr' - r'^dcj)' (g) d(j)' - dz' (g) dz' . (83) 

We immediately read from Eg . ((821) that 

5oo = (1 - v^) = l/7^ g'w = -^^'i = ^2/': 520 = = ^a;', g'^Q = 0. (84) 

After this (long) digression we return to Eg.([57| and Eg. (|60l) that in coor- 
dinates (x'^) which is (nacsjV), can be immediately written in the engineering 
format (of vector calculus) as: 

^**Pay attention with the notation used. 
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D' = eE', 



Then, we finally get 



D' = eE', 



1 



(85) 
(86) 

(87) 
(88) 



The polarization and magnetization vector fields in engineering notation are, 
respectively, then: 



P' = {e- 1)E', 

M' = (1 - -)B' + -7^(1 - efi)v X E'. 
11 H 



(89) 
(90) 



To obtain the expression of those fields in the laboratory (in coordinates (x'') 
naturally adapted to i = eo = d/dx'^ = d/dt) it is only necessary to recall that 



, 5x« dn!^ , dx" dx^ , ax" 9x^ 

F„„ = TT^TT^Fa/j, H = 7^377 ^ZT;: Ha/3, 



Now, from Ea. (|75|) we have 



na,3. (91) 



/9x"\ 



Then, we have, e.g., 
K 

K 
K 

B' 



1 

— wy cos Lot' 

u)x sinwi' 
\ 





— sinwt' 
cos Lot' 







1 / 



(92) 



Ex cos Lot + Ey sin ujt + lox'Bz 



—Ex sin Lot + Ey cos Lot 



WJ/'Bz 



= Ez - cjyBy - cjxBx, 
— Bx cos Lot + By sin Lot, 

— Bx sin w< + By coswi, 

B,. 



(93) 



We recall also that writing at time t in the laboratory the position vector 
of a point p = (x, y, z) in engineering notation as x = xCx + yey + zez and the 
(3-dimensional) angular velocity field of the frame V as u = fJe^ we have 



V = Q,xx. 



(94) 
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Finally we write the relation between the charge and current densities as ob- 
served in the laboratory and the rotating frame. From J — J'^d'^ — J^^^, with 
{J\J\J\J'') = and J'3) = {p\y^,j'y^) we have 



which gives 



J'^ = ^J", (95) 



p' = (96) 

= ojy'p + jx coswt + jy sinwi, (97) 

= — wx'p — jx sinwt + jy coswt, (98) 

j'.=h- (99) 



4 Jump Conditions for the Fields F and G at 
the Boundary of a Moving NDHILM 

The interface between a moving NDHILM described by the velocity field V 
and the vacuum defines a hypersurface E = in Minkowski spacetime. The 
jump conditions can be deduced from Maxwell equations and are well known. 
A very simple deduction of the discontinuity of the fields B, D and H is 
given in [10]. Here we write the jump conditions in differential form style for 
the case where the free current ★J = 0. Denoting as usual by [F], [G] and [★G] 
the respective discontinuities of F, G and ★G at the boundary of medium and 
vacuum we have: 

[F] A dSls^o = 0, (100) 
[G]LdS|3^o = 0. (101) 

Observe that [G]L(iS|-^Q = implies, of course 

*[[G]LdS|3^J = - * [[G]^dS|s^J = [*G] A dSls^o = 0. 

Also, [F] A dSls^o = A [i^lls^o — ★ F)]^^^ — and we can write the 

jump conditions also as: 

[*F]^dE\^^, = 0, (102) 
[★G] A dSl^^o = 0. (103) 

Now, in the (nacs|L = ep) (x^) we have 

dE 

n := dE = — = Ua'd'' = 7i^i9„. (104) 

Definj^ n :— — the spatial vector field which is normal to the moving 

boundary S(t,x,y, z) = 0. Now, each spacial point of the moving boundary at 

^^The minus sign is necessary due to the signature of the metric. 
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time t has Newtonian velocity v. Observe that during a time interval Ax" = At 
a point at the boundary with arbitrary coordinates (x, y, z) will arrive at the 
point (x + WxAt, y + VjAt, z + At). The hypersurface S at (t + At, x + 
WxAt, y + WyAt, z + VzAt) will satisfy 

S(t + At, X + WxAt, y + VjAt, z + v^At) 
dE 



"(t,x,y,z) 



Ax^ = 



At 



dE 
dt 



(t,x,y,z) 
dE 



VxAt 



(t,x,y,z) 



(t,x,y,z) 



A c>E 



zAt 



dE 



(t,x,y,z) 



= 0. 

(t,x,y,z) 

(105) 



Thus, we get in engineering notation taking into account that n — — VS that 
Eq. (|105p implies 

^-nmv = 0. (106) 
at 

Then, we can write the jump conditions in its usual engineering notation as 

n»[S]=0, {n»v)[B]+nx[E]=0 (107) 

(108) 



r I I 

n»[D]={), {n»v)[D]+nx[H]=0. 



I I 

Exercise 11 Show that Ea. (jlOip implies the formulas in Eq. (|108p . 

Solution 12 We calculate [G\^dE\^^Q ~ [G]^n\^^^ = '^-'[G'll^^o = ^- ^ow, 
from Eq.(|5ip we have 



njG = nj(Z A £)) - A H)] 



and since l_sD — ^^_s{Dij'^) = we get 



nj(l A r>) = {njl)D - n A {IjD) = 71qD 



which in engineering notation reads 



Also, 



(ft • v)D. 

nj[*(Z A H)] = ir{n A I A H) 

= *{n,Hj-f' A 7" A 7^') 
= -n,;i7j*(7° A7* A7^') 



(109) 



(110) 



(111) 



(112) 
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which in engineering notation read: 

-nxH. (113) 

Using Eq. (|llip and Ea. (lll3p permit us to write the equation rij[G]\~^Q — in 
vector calculus notation as 

{n»v)[D] + nx [H] ==0. 

To show that nj[G]\~^Q = implies also n» [D] = it is enough to observe that 

Zj(nj[G]) = {lAn)4G] 

= n,(70Ay)j[iG"^„A7/j] 
= -n.(7"Af).[iG^^A7/j] 

2 ' V 7 • 7a 7 ■ 7/3 / 

-n,[G°']=0. (114) 

But G"* = -D' and thus we can write lj{nj[G]) = in vector calculus notation 
as nu [D] ~ 0. 

5 Solution of Maxwell Equations for the Wilson 
&; Wilson Experiment 

We now show how to find the solution of Maxwell equations 

dF = 0, d*G'=-J (115) 

for the famous Wilson & Wilson experiment of 1913. Here that experiment is 
modelled as follows: a cylindrical magnetic insulator of internal and external 
radii r^ and r2 respectively and which has uniform and isotropic electric and 
magnetic permeabilities e and /i is supposed to rotate with constant angular 
velocity lo in the z direction of an inertial laboratory i = Cq where there exists 
a uniform magnetic field Fq — ^^drhdc^ (or in engineering notation Bq = BoGz). 
Figure 1 illustrate the situation just described. 

We start by introducing some useful notation. First we write the Minkowski 
metric g as 

g^n^.uO'^^e'' (ii6) 

with 

6)° = dt, = dr, = rd(j), = dz. (117) 
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Figure 1: Details of the Wilson & Wilson Experiment. A magnetic insulator 
rotates with uniform angular velocity lj = ujEz in an external magnetic field 
Bo = Bo&z. An electric field Ei — /(r, Bq, w, e, /x)er is observed. 
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For our problem the moving boundary of our material has equations 5i = 
r — ri = and S2 = r — r2 = 0, so the normal to this surface is (E = 6^. For 
our problem we have J = 0. To proceed we write as an ansatz the solution for 
the electromagnetic field in the interior of the material 

Fi=£e" AO^ +B9^ AO^ (118) 

where £ and B are supposed to be functions only on the coordinate r. 

We recall from Eg . ([75)) that the 1-form physically equivalent to the velocity 
field V is 

V = 7(6>° + vO^) 

with V — uir and 7 = (1 — v^)^^/^. We now must solve the equation d* G = 
with the boundary conditions given by Eqs. (ll02p and (|103p . To calculate Gi we 
use Eg. ([5^. First we calculate v_iFi. We have 

^j{£ + vB)e\ 

and 

V A (vjFi) = 7^(f + vB)e°^ - 7\(f + vB)e^^. 



Thei 
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_ g(/i£-v^) + gv(M£- 1) ^01 ^ £Al-^le) + B{l-^iev') ^,^ 
m(1-v2) + m(1-v2) 



Since 



*{e° A 0^) = (6»" A 0^)jTg = -6/2 A 
A 0^) = (0^ A 6>2)jTg = 0° A 0^ 



we have 



* Gi = -O^^ + £0°^ (120) 

Then d-kG — gives 

_^^123 _ ^gl23 _ ^^013 _ Q 

and so our problem resumes in solving the following two trivial ordinary differ- 
ential equations 

die K- ^ dC , . 

— + - = and — = 0, 121) 
dr r dr 

with solutions 

/C = — and £ = C2 (122) 

r 



^It seems to have a misprint sign in the formula appearing in [T]. 
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respectively, where ci and ci are integration constants. So, we have 

r 

r 

Now, if we recall that outside the magnetic insulator we have -kGo = *Fo = 
Bq^"'^ we get using the jump condition Ea. (ll03p that 

(★Go - *Gi) AO^ ^0, (123) 

i.e., 

- Bo0°" + C20°" - -e'^^ = (124) 
r 

from where it follows that the functions /C and C are 

/C = 0, C = Bo. (125) 

Then from Ea. (|119p we have the following system of linear equations for the 
functions £ and B: 

r £{fie-v^) + Bv{fie-1) = 0, 
I Ki-^'') ^ ° 



whose solution is: 



^-^o^^^ff^y B^-^J-t^f^^. (127) 
e(l — w^r^) e(l — cj^j'^) 



So, finally we have 



a;r(l - /ie) (w^r^ - ^e) 



e(l — w^r-^) eyl — uj^r'^) 

and the electric and magnetic fields and ifi as determined in the labo- 

i I 

ratory frame inside the material are 

^. = ^°.F. = Bo "f -f,> \ (129) 
I e(l — w^r^) 

ffi = 0". ★ = -Bo^^^J^— (130) 

From Eq. (1129^ it follows immediately that the potential V between the internal 
and external parts of the material as shown in Figure 1 is, when w^r^ << 1, 

y = Bo^i(i-M£)/;>rf^ 

= ^B,oj-^{l - fie){rl - rl), (131) 
which is the value found by Wilson & Wilson [21] . 

Exercise 13 Calculate P , M , the bound current J and the bound charge p. 

II I I 
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6 Extracting Energy from the Magnetic Field 



We now propose a way to extract energy from a magnetic field using the results 
just obtained above. In order to do that we first recall that when the magnetic 
insulator in Figure 1 which has momentum of inertia / is put in rotation with 
constant angular velocity uj = cje^ it acquires an angular momentum 

Lmec = I^- (132) 

For the preliminaries theoretical considerations in this section we suppose that 
there are no energy losses due to friction of the rotating magnetic insulator with 
its supporters nor losses due to Joule effect on electric wires. As a consequence 
the total angular momentum of the system (i.e., the rotating dielectric plus 
the electromagnetic field) must be conserved. We suppose moreover that the 
system does not dissipate energy through radiation. Under these conditions the 
total angular momentum of the system is 

(133) 

where by [12] Abraham's formulsF^ 

Leiec = Jdxdydzx x (Ei x Hi). (134) 

We choose our coordinate system such that the origin lives in the middle of the 
rotating axis of the dielectric. Under these conditions 

E:xH: = _B^ Ml-M.)^K^^r^ -Me) .^ ^ „^ 
e^(l — Lu^r^)'' 

Then 

X X (£;i X Hi) = (ze^ + rer) x f{r, Bo, /i, e)e^ 

= z/(r, Bo, w, ^, e)ez x + r/(r. Bo, w, ^, e)er x 
= -z/(r. Bo, w, e)er + rf{r, Bo, u>, n, e)e^ 
= -z/(r, Bo,w,Ai,e)(cos0ex + sin^Cy) + r/(r. Bo, w, ^, e)ez. 
Thus the integral corresponding to the component vanishes and we have 
Leiec = r^/(r. Bo, w, A*, e)drd4)d-z 

= 27rZe2 J^^'r^f{r, Bo, n, e)dr. (135) 
where Z is the height of the cylindrical dielectric. When w^r^ << 1 

= -^B^c.(^£-l)(r^-rf)e.. (136) 



^^At first sight it may seems strange that a magnetic and an electric field coming from 
different sources may storage angular momentum. However that they do is an experimental 
fact as showed for the first time by Graham and Lahoz only in 1979 [4l[7]. 
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It is a remarkable fact that Lf-iec points in the opposite direction of imec if 
/ze >1 (as it is the case in the Wilson & Wilson experiment according to [TT] . 
Since the total angular momentum of our system must be conserved we see that 

Lmec must increase if there is no energy losses due to friction and Joule effect 

in the wires and we could use the potential V given by Eq. (|13ip to power an 
electric machine. Of course the energy powering the electric machine can only 
be coming from the energy stocked in the external magnetic field. It will work 
while the external magnetic field is presented. Of course extracting energy from 
the magnetic field will make the value of Bq to decrease and with this decrease 
the potential V will also decrease and at the end the machine will stop working. 

Remark 14 If <I the electromagnetic angular momentum is in the same di- 
rection of the mechanical angular momentum, which must start to decrease after 
the device is put in rotation. Then, the electric field will decrease also making 
the electromagnetic momentum to decrease whereas the mechanical angular mo- 
mentum must then start to increase again. It is not clear at this moment if 
the system will oscillate between a minimum and maximal mechanical angular 
momentum or will stabilize. 

Remark 15 Real World Machine. In the real world where friction and Joule 
effect are always present in order to extract energy from the magnetic field we 
need a machine like the one in Figure 2 where the potential difference V given 
by Eq. (jl31l) is used first, if necessary, to power in a compensator motor just 
enough to restore eventual losses due to friction and Joule effect on the electric 
wires, and secondly V is also used to power a motor to generate useful work. 



7 Conclusions 

In this paper we recalled how to correctly solve Maxwell equations in order to 
find the electric field in the famous Wilson & Wilson experiment using the theory 
of differential forms with permits an intrinsic formulation of the problem. We 
present also the correct jump conditions for the F and G fields in an invariant 
way when the boundary separating two media is in motion. We give enough 
details for the paper to be useful for students and researchers. Moreover we 
use the theoretical results to present a surprising result which for the best of 
our knowledge is new: a machine that can extract energy from an external 
magnetic field. Under appropriate conditions (/ie >I) the machine once puts 
in motion will start increasing the electromagnetic angular momentum stocked 
in the electromagnetic field in the direction contrary to the original mechanical 
angular momentum of the device which thus will start increasing (theoretically) 
its angular speed until Lor2 = 1, thus generating a big potential even in a small 
external magnetic field which can be use to produce useful work. 

Of course, a machine like this one if in orbit around a neutron star can 
produce a lot of energy in a simple way than the hypothetical machine described, 
e.g., in page 908 [9] projected to extract energy from black holes by throwing 
garbage on it. 
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A Some Useful Formulas 



(138) 



In this appendix we present some useful identities of the theory of differential 
forms with has been used several times in the main text. 

For a e sec l\^T*M and Ar € sec /\''T* M, Br € sec /\''T*M it is 

a^Ar A Br) = (ajA,.) A Br + Ar A {a_iBr). (137) 
Ar A i^Bs — Bs A -kAr; r — s 

9 9 

Ar ■ *-Bs = Bs ■ ^Ar ; r -\- s — n 

9 9 

Ar Ai.Bs = {-ly^"^^^ -k{ArABs)\ r<s 

9 9 

ArJ-kBs = {-lY'-k{ArABs); r + s<n 
a 9 

'^j\.f — y4.j-'_iTV7 
9 " 

*Tg = sgndetg; ^1 = r„. 
9 g 

Acknowledgement Authors are grateful to Flavia Tonelli for the figures. 
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Figure 2: Machine to Extract Energy from an External Magnetic Field. 
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